Microscopic Oscillations in the Quantum Nucleation of Vortices 
Subject to Periodic Pinning Potential in a Thin Superconductor 
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We present a theory for the decay of a supercurrent through nucleation of vortex-antivortex pairs in 
a two-dimensional superconductor in the presence of dissipation and of a periodic pinning potential. 
Through a powerful quantum electrodynamics formulation of the problem we show that the nucle- 
ation rate develops oscillations in its current-density dependence which are connected to the pinning 
periodicity. A remnant of the dissipation-driven localization transition is present, and an estimate 
of the nucleation rate suggests that these effects might be observable in real thin superconductors. 
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There has been a strong revival of interest, recently, 
in the physics of magnetic vortices in type II and high- 
temperature superconductors Phenomena involving 
such vortices pose well-defined problems which are of 
great importance in the technological applications of the 
new high- Tc materials. Of particular interest are the stat- 
ics (the vortex-glass transition) and dynamics (flux creep, 
vortex tunneling) of vortex systems in the presence of 
pinning impurity centers and an externally applied su- 
percurrent. Most of these issues are connected with the 
residual resistence due to the thermal or quantum mo- 
tion of the vortices and as hindered by the presence of 
pinning centers. 

In a recent article, we have pointed out that such 
residual resistence should develop even in the absence of 
an applied magnetic field. This is due to the sponta- 
neous nucleation in the bulk of vortex-antivortex pairs 
induced, in the presence of a supercurrent driven by 
an external source, by the electromagnetic-like Mag- 
nus force acting on the quantum-mechanical vortices 
thought of as electron-positron- like pairs. This quantum 
nucleation occurs homogeneously throughout the sam- 
ple. Exploiting this analogy with quantum electrody- 
namics (QED), a powerful "relativistic" quantum field 
theory (QFT) approach has been set up to study vor- 
tex nucleation in the two-dimensional (2-D) geometry 
and in the presence of quantum dissipation and of in- 
finite barrier-height pinning centers. The central re- 
sult is that the nucleation rate P has a strong exponen- 
tial dependence on the number current density J, given 
by PocTyi/^r^eZ/J"^ exp{-77e//£:^^/47rj2}. Here, rjeff is 
an effective viscosity coefficient as renormalised by the 
magnetic-like part of the Magnus force, and Eor is the 
rest or nucleation energy of a single vortex as renor- 



malized by screened Coulomb interactions and (fake) 
Landau-level corrections. This exponential dependence 
would make the vortex nucleation (folded, e.g., into the 
sample's resistance) observable in a rather narrow range 
of J- values (see below). We stress that whilst the possi- 
bility of homogeneous vortex nucleation as a decay mech- 
anism for a supercurrent was independently and first pro- 
posed by Ping Ao, |3| the extension to dissipative dy- 
namics, leading to qualitatively different results, was first 
treated in our recent paper, to which the interested 
reader is referred for details of the formalism. 

In this Letter we extend our formulation to the im- 
portant situation in which finite barrier-height pinning 
centers are present and we initially envisage a geome- 
try in which the pinning potential is periodic, with a 
lattice spacing d = 27r/fc and finite amplitude Uq. We 
estimate the order of magnitue of the production rate P 
to find that under reasonable conditions the nucleation 
phenomenon ought to be observable. Moreover, we show 
that a remnant of the dissipation-induced vortex localiza- 
tion transition |^,|| is present through a J- dependence of 
P in which microscopic (or, depending on d, mesoscopic) 
oscillations are noted, due to the underlying periodicity 
of the pinning potential, and might also be observable 
under suitable conditions. Beside spontaneous bulk nu- 
cleation in the absence of a magnetic field, our formalism 
and results can also be applied to vortex tunneling from 
the edge, extending the work of Ref. |^-|^ with the aim 
of obtaining an explicit formula, complete of prefactors, 
for the tunneling rate. Details of our calculations for the 
periodic pinning potential and the treatment for the edge 
tunneling problem will be presented elsewhere. 

We begin by summarising the problem at hand, essen- 
tially a single- particle one, through the classical equation 
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of motion for the vortex moving in the (x — y) plane of 
a 2-D superconducting film with a current flowing in the 
y-direction: 



mq : 



-VU (q) + eE • q — eqxB - r/q 



(1) 



Here, m— >0 is the inertial mass of the vortex of topolog- 
ical charge e = ±1 and 2-D coordinate q(i) subject to 
the Magnus force having an electric-like E — 27rJx as 
well as a magnetic-like B = 27r^psSZ {ps being the su- 
perfluid 3-D number density and s the film thickness). 
The e.m. current density J*^™ flowing through the sam- 
ple's cross section is related to the number current den- 
sity via J = J*^™s/2e. We work in the units system in 
which Ti = 1. Also, is the friction coefficient and U{q) 
the phenomenological periodic pinning potential which 
for convenience we take to be of the cosine form with 
wavenumber Q = (/c, fc): 



C/(q(i)) =C/o^[l-cos(fcg,(t))] 



(2) 



a=l 



The quantum-mechanical version of Eq. (|I]) is attained 
through the Feynman path-integral (FPI) transposition 
in which dissipation is treated quantistically through the 
formulation due to Caldeira and Leggett. |^ Particle 
nucleation (closely related to tunneling) is dealt with 
through a scalar-QED, 2-D version of Schwinger's calcu- 
lation for the vacuum decay rate T in the presence of the 
crossed uniform fields (E,B). The formula for T reads 



r _ 2 



Im 



dr 



Vq{t)exp{- I dtCs} 



(3) 



where L^T, with D = 2 here, is the space-time volume 
of the sample, Sq is the vortex nucleation energy, and the 
Euclidean single-particle Lagrangian for closed trajecto- 
ries in space-time is 
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Ce = -jm^cif^q^, - -iq^Ff,,,q^ + F(q) 
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mk±l + ^rukujl ( Xfc + — ^q 
2 2 \ mkUJl 



(4) 



Here, the set of {x^} represents the Caldeira-Leggett 
bath of harmonic oscillators simulating quantum friction. 



-5{uj — ujk) = rjLo ex'p{—co/n) (fl being 



with f yit ■ 

some frequency cutoff) fixing the spectral distribution 
of frequencies (we restrict the present discussion to the 
ohmic case). Also, F^^, is the uniform Magnus field tensor 
and V^(q) = 2£[)U (q) is the relativistic QFT counterpart 
of the classical potential, Eq. (||), thus having amplitude 
Aq = 2£qUo (we have assumed Aq/£q<^1). Also, 7713 = ^ 
and mi = m2 = ^ with ^ — £0/ m^oo playing the role of 



the square of the speed of light in this formalism, where 
the non-relativistic limit corresponds to the overdamped 
case (which we consider). 

Before delving into the evaluation of the FPI in Eq. 
(^), we briefly dwell on the physics of the problem at 
hand in the absence of the Magnus force, and on the 
meaning of Eq. As is well-known, the motion 
of a quantum particle in a periodic potential in the pres- 
ence of friction displays a localization transition whereby 
for rj < rjc the particle is mobile whilst for rj > rjc it is 
confined, rjc = fc^/27r = 2-K/d^ being a sharp threshold 
independent of the potential's amplitude Aq. This is re- 
produced by our formalism, with E = B = and 7^0 in 
Eq. (^ , with the confined phase treated within the varia- 
tional "blocked-renormalization" scheme first introduced 
by Fisher and Zwerger, through which a "mass" is 
dynamically generated for the damped modes' propaga- 
tor. When the Magnus force-field is turned on (at least 
its electric-like part) a pole develops in the bare propa- 
gator, signalling the onset of mobility for all values of 77. 
This divergence is in turn embedded in the formula for 
the pair production rate, Eq. (^. This one comes from 
a formulation of vortex nucleation |^ in terms of a scalar 
quantum field 4> with QFT Lagrangian 

£(0) = Da<jy*Do(j) - \t>^* ■ - {£l + y(q)) 0*0 (5) 

where 13^ = — ieA^ with = ^F^,^q^, in which the 
process of particle-antiparticlc pair production is quite 
naturally encoded. In this QFT formulation, the vacuum 
decay rate is computed through the FPI on all closed 
loops in space-time weighted with the single-particle Eu- 
clidean Lagrangian corresponding to (^). 

We now come to the evaluation of P, in the presence of 
the periodic potential which calls for a renormalization 
group (RG) type approach. Indeed, apart from the sub- 
tle question of dynamic mass generation for rj > rjc, the 
problem at hand is that of a sine-Gordon-like model in 
(1-1-0) dimensions. Integrating out the Euclidean "time"- 
like component (73 (i), we reach a formulation in which the 
electric-like and the magnetic-like Magnus field compo- 
nents are disentangled. In terms of Fourier components 
q„ = q(t^„), with lj„ = 2Tm/T: 

rT ^ ^ ^ 

/ dtCsiq^,)^ dtCsici) = T^^i-^fujl + -rj\uj^\ 

-E^Sai}qai(^n)qai-(^n) + / dtV (q) + magn. part (6) 

Jo 



In this paper, we take the stance that the only role of 
the magnetic-like field in the presence of dissipation is 
to renormalize the nucleation energy, fo~^^oi? 



- ^0 



^ rf -\- B'^ , as well as the friction coefficient, r/- 
{rf' + B'^)/ri (denoted by f] in what follows). This was in- 
deed one of the main results of our investigation, ||^ and 
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it amounts to an effective reduction of the D-dimensional 
problem to a one-dimensional system in the presence of 
the sole electric-like field. [|| A crucial role is played by 
dissipation, of course, as otherwise we would be con- 
fronted with the Azbel-Hofstadter problem in the pres- 
ence of both periodic pinning and magnetic field. The 
evaluation of the FPI proceeds by means of the inte- 
gration over the zero-mode, qo, and the high-frequency 
modes, qn with n > 1, where we take the point of view 
that, as for the free particle in the uniform field, |Q the 
leading term in F comes from the divergence of the FPI 
associated with the lowest mode coupling to E. We apply 
this scheme to the formula for the rate, Eq. (^, where 
the effective 1-D Lagrangian is the same as Eq. (^) , with 
7 = and renormalized Sq and r/, and taking qi=q to be 
strictly a one-dimensional coordinate. 

The effect of the modes g„ with n > 1 is taken 
into account through a frequency-shell RG-like argu- 
ment by which, for the interesting rj > rjc case, these 
modes simply renormalize Aq and are integrated out 
as effectively "massive" [0, that is via a Lagrangian 
'C(gn) = T (i77|w„| -I- — E'^) q„q^ (we will see below 
that at the relevant singularity = ^r]\uJi\, so that ef- 
fectively n^n — 1). The result of this integration is an 
additional, entropy-like renormalization of the activation 
energy, to be included in Sqr, beside the renormaliza- 
tion of Aq. Notice that in the mobile phase, t] < r]c, 
the pinning potential's amplitude would renormalize to 
zero, suppressing the oscillations (to be discussed below) 
in the vortex nucleation rate. We are therefore left with 
the integration over the modes qo and qi , with a potential 
term 



dtVR{qQ,qi{t)) ^ Aqr / dt cos k{qQ + qi{t)) - Aqt 



iAQRTJo{2k\qi\) cos kqo - Aqi 



(7) 



In this equation, Jq is the Bessel function and the new 
amplitude Aqr, as renormalised by the "blocked" n > 1 
modes in the confined t] > rjc phase, is given by 



Aqr = 



A 



(p+rt--i) 



if ?7 < 7/c 

if J] > T]c 



(8) 



where we take the cutoff integer n*{T) — [ilr/27r]::^l and 

a dimensionless mass /i — = ( ^q^^" j has 

been introduced. Notice that, from the variational ap- 
proach, /X = for r/ < rjc- Also, in Eq. (^ the phase of 
the qi mode has been integrated out (a procedure gener- 
ating an expansion controlled by the non-divergent am- 
plitude Aqp), allowing us to integrate out the go-mode 
exactly. This leads to the expression 
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— ^=Im I dr7V(T)e"(^««+^'')^ x 



d|gipe-(-''-^'^)l^"^l'/o(^oflTJo(2fc|gi|)) 



(9) 



where M{t) is a suitable normalization factor and /q is 
the modified Bessel function. The vortex nucleation rate, 
therefore, bears some sign of the friction-driven localiza- 
tion transition, through Aqrt as given by Eq. (||). It is 
clear that the singularity in r from the gi -integral occurs 
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nrj/E ; taking a suitable ie prescription H we 



AA(ri) 



2L2 47rJ2 



V47rJ2 



(10) 



where we have set E = 2nJ . The presence of the 
Jo(2fc^Ar) argument in the correction factor due to the 
pinning lattice thus gives rise to oscillations in F (hence in 
the sample's resistence) through the parameter 2fc£jv — 
Airtisi/d- Vortex nucleation (or tunneling) is therefore 
sensitive to the corrugation of the pinning substrate. 

In order to compute the nucleation length we 
have used a saddle-point procedure in Eq. (^, notic- 
ing that for large r: Jo(xr)^ , el^l"^. This leads to 

y/2'ir\x\T 

a self-consistent equation, -|- — ^oK|^o(2A:£Ar)| — 
47r2 j2£2^, which we solve to first order, neglecting correc- 
tions for small J: 



1)2 ^ ^QR + ^0 



tjv- 



47r2j2 



47r2j2 
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OR 



An 



nJ 



(11) 



As for the normalization factor A/'(ti), we get (e is here 
Euler's number) 



-^(^i) = T^"*(^i)(l + + M«*(ti))(Vti)'/^' 

(12) 

Our final expression for the rate F (with In given by 
Eq.([l])) is 

F = FoX(J) 



^2 327r2j 



KiJ) 



(13) 



where K{J) is the oscillating correction factor due to the 
pinning substrate. 

So far, the computation referred to T = 0. In 
we proposed to take into account finite temperatures by 
adding incoherently, in the above formulas, a term pro- 
portional to T (note ks = 1): -I- ^^T, such that 
for J— >0 the standard Boltzmann exponential form in the 
vortex density, pyOcy/T, is recovered in the absence of pin- 
ning. This might be a reasonable qualitative treatment of 
thermal effects for temperatures less than KiSnJ^ /rjEpR. 
Furthermore, we propose, following Minnhagen ||l^, to 
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FIG. 1. Vortex-antivortex production rate (in iim~'^s^^, 
logj^Q scale) in the presence (F) and absence (Fo) of the 
periodic pinning potential, versus current density J*^™ (in 
A cm~^ X 10^). K{J) is the current-dependent correction 
factor due to the periodic substrate. 



duction rate, both in the absence and in the presence of 
the above-specified periodic pinning potential. We stress 
that, despite our possibly optimal choice of parameters, 
we believe that the order of magnitude of F is such as 
to make vortex nucleation observable within acceptable 
values for the current density of some high- Tc materi- 
als. When pinning is present, we predict for the first 
time the possibly observable oscillations of r( J), reflect- 
ing the underlying pinning periodicity. This should pro- 
duce microscopic oscillations in the measured resistence 
at relatively low temperatures (T = 5 K in Fig. 1). For 
higher temperatures, vortex nucleation would be essen- 
tially thermally activated and therefore outside the scope 
of the present investigation. 

We hope that these findings will stimulate further re- 
search, both theoretical and experimental, on the quan- 
tum mechanics of vortices in thin superconducting films. 

One of us (GJ) is grateful to the International School 
for Advanced Studies in Trieste, where part of this work 
was carried out, for hospitality and the use of its facilities. 



account for vortex-antivortex Coulomb-like interactions 
by means of a current-dependent (and temperature- in- 
dependent) activation energy: £r{J) = £{)R.\Q-{Jmax / J), 
with £oi?, including all other renormalization effects. This 
expression entails an £r infinite for J = (as is appro- 
priate below the Kosterlitz-Thouless (KT) transition in 
the film) and changing sign, thus becoming unphysical, 

Jjnax • 

For illustrative purposes, in order to estimate the pa- 
rameters involved in our final expression for the rate 
F (Eq. (p^), we take Jmax an order of magnitude 
higher than the single-crystal YBCO critical current, p^ ] 
namely J^™^ = 10^ A cm'^, and take £flw80 K for 
J = 10^ A cm^^, of the order of magnitude of the KT 
transition temperature [Q. The important friction coef- 
ficient is taken to be ?7wlO~^ A , from the Bardeen- 
Stephen formula applied to YBCO films. A most 
sensitive parameter is related to the amplitude of the 
pinning potential, e = Aq/S^j^, which, like the pinning 
lattice spacing d, is entirely unknown. We have taken 
a negative e = —0.5, borrowing from classical nucle- 
ation the point of view that vortex production is actually 
aided by the presence of pinning centers. Finally, c? = 50 
A concludes our illustrative case which is clearly of the 
confined-phase type. As > r]c, we have taken and 
AoRfnAo, for simplicity's sake, and find no significant de- 
pendence on the most uncertain of all parameters, Q. Its 
dimensions being the same as £q, in Fig. 1 we have taken 
n = fpfl/e. Also, we take the film thickness s = 10 A as 
the typical interlayer spacing. This figure reports the pro- 
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